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A specific scale factor in Robertson- Walker metric with the prospect of giving the 

overall cosmic history in a unified picture roughly is considered. The corresponding 

' energy-momentum tensor is identified as that of two scalar fields where one plays the 

(N ■ 

roles of both inflaton and dark matter while the other accounts for dark energy. A 
preliminary phenomenological analysis gives an order of magnitude agreement with 
i observational data. The resulting picture may be considered as a first step towards 

^ ^ a single model for all epochs of cosmic evolution. 
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I. INTRODUCTION 



There is an intense on-going research to understand the natures of late-time acceleration 
1| (whose standard explanation is dark-energy j2|), dark-matter jsl, and the inflationary era 
^, |5|. A detailed and deflnite formulation of each of these issues by its own is essential and 
very important for the future direction of cosmology. However how to relate these in a single 
formulation and unify ally eras of cosmology (namely, inflationary, radiation dominated, 
matter dominated, the current late-time acceleration) is as essential as the study of each 
era separately. This is not only due to the fact that we must eventually put all these into 
single picture but it is necessary for a better and correct formulation each of these issues. 
This paper is an attempt in this direction i.e. to obtain an overall picture of cosmic history 
in a single model. Some of the other studies in this direction may be found in Q-lsl. I hope 
the study given here is simpler and more concrete while being minimal and formulated in 
standard framework i.e. two standard scalar flelds in the usual 4-dimensional Robertson 
Walker metric and the usual Einstein-Hilbert gravity. 

In this paper I consider a speciflc scale factor in the usual 4-dimensional Robertson- 
Walker metric. The scale factor is chosen in a such a way that it has a prospect to account 
for inflationary, matter dominated, and late-time acceleration eras of the universe. Then I 
check this expectation. First I flnd the corresponding energy-momentum tensor and identify 
it by that of two scalar flelds. The flrst one mimics inflation at very small times and then 
mimics (dark) matter at intermediate times. The second one is identifled by dark energy. 
Hence this model accounts for all epochs of the universe except the radiation dominated 
one. The content of this universe is similar to our own except it does not contain baryonic 
matter and radiation. This universe is similar to our own, given the fact that the present 
ratio of baryonic matter density and radiation to the total energy density of the universe 
is 4 % and hence negligible and remains negligible (at gravitational level) in most of the 
cosmological evolution except in the radiation dominated era. Then I use cosmological data 
to constraint the parameters of the model and apply these to some redshifts and to the 
corresponding time data to check the phenomenological viability of the model. There is an 
order of magnitude agreement with data. In my opinion the results are encouraging to look 
for a more elaborate form of the model where baryonic matter and radiation are included, 
and where a more thorough study of the parameter space is investigated. 
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II. THE MODEL 

We consider the Robert son- Walker metric 

ds^ = g^ydx^ dx" = —dt^ + a^{t)gijdx^dx^ (1) 

We take the 3- dimensional space be flat, i.e. Qij = 6ij for the sake of simplicity, which is an 
assumption consistent with cosmological observations [1] . The key assumption in this paper 
is the following ansatz 

a{t) = [pi + p2a2t] Exp[-bi{a2t)-^^^] (2) 

where pi, P2, cl2, bi are some constants that to be fixed or bounded by consistency arguments 
or cosmological observations. The corresponding Hubble constant and its rate of change are 
given by 

7 

^ _ a _ a2[6p2(Q2^)6 + + P2a2t)] 

« Q{a2t)^ (pi + p2a2t) 

il_ ^ (a2t)^[36p^(a2t)^ + 7&i(pi + p2Q2t)^] 
H"^ [6^2(02^)6 + h{Pi+P2a2t)f 

and the acceleration of the scale factor is 

a _ hi[hi{pi + P2a2t) + (a2^)^(-7pi + 5p2Q2t)] 
a 36t2(a2t)3(pi +^20-2^^) 

where the dots on top of the letters stand for time derivative. 

The following observations about the scale factor a{t) are in order; One notices that a{t) 
is positive for all values of t provided that 

Pi,P2,a2,bi > (6) 

^ is positive for extremely small values of 02^, where the leading term in ^ is the pi term; 
^ is negative for the intermediate values of a2t, where the leading term in the numerator is 
the — 7(a2t)^ term; and ^ is positive again for the larger values of 02^. Note that the present 
era corresponds to very large values of t, not the infinite value of t where the acceleration 
is zero. One may see the general form of the evolution of ^ for a set of phenomenologically 
relevant parameters in the next section in Figure [2j Moreover it is evident from Eq.(|l]) 
that here is almost zero (i.e. slow- condition is satisfied) if a2t is taken sufficiently small. 
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Therefore the scale factor ansatz given above, at least in principle, is suitable to account 
for all four eras of cosmic expansion; inflation, radiation dominated era, matter dominated 
era, and current accelerated expansion era. In the analysis given below first I will determine 
the Einstein tensor and the corresponding energy-momentum tensor. I will identify this 
energy-momentum tensor with that of two scalar fields. Then, after using phenomenological 
considerations, the parameters (i.e. pi, p2, 02, &i) are numerically constrained. I will check 
the phenomenological viability of the model. It will be seen that the scalar fields may be 
identified by inflaton, dark energy and dark matter, and the corresponding picture is that 
of a universe that consists of only dark energy and dark matter (that also serves as inflaton 
at early times). Given the fact that, at present, more than 96 % of the universe consist of 
dark energy and dark matter this universe will be considered as a universe that is similar to 
our own in its overall cosmic history except in the radiation dominated era. Although the 
results obtained here have only order of magnitude agreement with observations, the results 
are encouraging for adopting this model as a starting point for a more elaborate formulation. 

The components of the Einstein tensor for the metric given by ([1]) with the scale factor 
in Eq.(|2]) are 



Goo = = ^^2[6p2(a2t)^ + hjpi + P2a2t)]y 
12(02^)3(^1 +P2a2ty 

i^^ , Tji\ 36p2(a2t)i + 362(pi+p2a2t)2 + 26i(a2t)^(-7p? + 4pip2a2t + llpiap') 
Gi,- = -(2- + ii = — g^^ 

« 36t^(a2t) 3 (pi + p2a2ty 

(8) 



Provided that we identify the source of the energy-momentum tensor as a collection of n 
real scalar flelds, its general form is 



n 1 " 

T^,u = ^[<9^0i<9j.0i + gf,u[--g^^ ^ dr(i)idp(f)i - V{(f)i, 02, ■, 0n)] 

i=l i=l 
n 

Too = H = Y, 7;^l + Vi<Pi, 02, -, 0n) (9) 



2 

i=l 



1 • 



Tij = [ 5Z o*^' - '^^^ 0n) ]gij (10) 



2 

i=l 
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After using the Einstein equations, we make the identification 

Gil o r^trr , 36p|(a2t)^ + 76i(pi +^2^2^)^ 



Gy^l = G,, + ^ = 8nG{Too + ^) = ""^^^"^^^ T ^^^^ " ^ ^ (11) 
^ ^711 9n m\a2t)l{pi+p2a2tY 



i=l 



IQtiGV = GQa - ^ = STrGfToo-— ) 
911 911 



18t2(a2t)^(pi+P2a2t)2 



(Pi+P2a2t)2 6{a2t)i 18(a2t)^ (a2t)^(pi + ^2^2^) 
Eq. ( ITT]) may be used to identify the scalars that act as the source of the Einstein equations 

Mt) = ci(a2t)-^ (13) 
(f)2{t) = C2 In {pi + p2a2t) (14) 

Writing the potential V in terms of these fields and satisfying the field equations 

dV ■ ■■ dV 

V^V^0i - 77— = -3i^0i - 01 - 77— = (16) 

901 (901 

dV ■ ■■ dV 

V^V^02 - 77— = -3i/02 - 02 - 77— = 17 

902 a 02 

identifies V as 

o /^T/ o 2 2 / 0*^2% , 1 ,2 2/*^l\28 2j, /'^1\26 

SttG^ = 2^202 exp -2— + —b^a^{—) - —a^bii—) 

C2 12 Ci 3o Ci 

+ P2a^6i(^)i^exp(-^) (18) 

Cl C2 

Then Eqs. (116111 71) are trivially satisfied for all values of the parameters, 02, 61, pi, P2- 

Next we will constrain these free parameters by phenomenological considerations and see 
if it gives a consistent and viable picture of the main lines of the cosmic history (except the 
baryonic matter and radiation). However before a phenomenological analysis it is necessary 
to identify which term in the above analysis corresponds to infiaton, which one to dark 
matter, and which one to dark energy. Before beginning the discussion it is worthwhile to 
note that both of 0i and 02 survive during all epochs of cosmic history. However only one 
of them is dominant in a given era of cosmic evolution. The infiaton term must be the one 
that is dominant and causes a huge cosmic acceleration at the time of inflation (i.e. at very 
small times). After examination of Eq.(I5]) and Eqs. (lll|12p we see that the dominant terms 
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(of huge contributions) for early times are proportional to (a2^)~^ oc (f)'l^. At intermediate 
times the dominant term is the term proportional to — 7pi in ([5]) i.e. the (02^)"^ oc term 
in ( lll|[T^ . Therefore 0i accounts for both of the inflationary and dark matter dominated 
eras. At late times the dominant contribution is due to the terms of the form 7 — — rw i.e. 
the terms containing 02. Hence 02 may be identified by dark energy. Because 0i is identified 
by dark matter its coupling to standard model particles must be small. However it must 
have large enough coupling with standard model particles to generate enough reheating. 
This may be accomplished by assuming 0i be electrically neutral and be color singlet. Even 
it may be taken to be a singlet under the whole SU{3)c ® SU{2)l ® U{1)l group of the 
standard model and couple to standard model particles indirectly (say through Higgs field) 
as in jo], Q. Another point would be a detailed study of the potential in (ITS]) especially 
to determine the effective range of 0i, in connection with its identification as dark matter, 
that is quite difficult due to the highly non-linear form of the potential. In fact all these 
points will arise when baryonic matter is included into the model and a more comprehensive 
and elaborate extension of this study is done in future. After these remarks we return to 
our main objective in the following paragraphs to check the phenomenological viability of 
the model. I put rough constraints on some of the free parameters, 02, &i, Pi, P2 through 
a rough empirical analysis. The cosmological eras that I employ to put constraints are 
the inflationary era, the present day, the onset of matter dominated era, and the time of 
reionization. I also consider the time of matter - radiation decoupling time. 



III. COMPATIBILITY WITH OBSERVATIONS 



The value of the present value of scale factor is taken to be one by convention. This 
implies 

1 = ao = a(to) = (Pi +^202^0) exp [-61(02^0)"^] (19) 

=^ (pi+P2a2to) = exp[6i(a2to)"^)] = /? > 1 (20) 

where /3 is some constant to be determined from observational data. We exclude the case 
/3 = 1 since it corresponds to infinite time for the present age of the universe. Next consider 
the observational values of the present value of the Hubble constant H{t = to) = Hq and 
the age of the universe t^. The observational values of Hq = , -,77^2 r 1 — TT^sHTr^' 



to = 13.69 ± 0.13 Gyr given by Particle Data Group (PDG) Q gives 

0.998 < Hoto < 1.018 (21) 

This implies HqIq ~ 1. Although the Hq and to values in f l2T|) are the most standard values, 
there are different observational values for Hq and to as well. For example Reese et. al. finds 



a value of Hq smaller than the PDG value by approximately 16% ll| although it may be 
ascribed to underestimation of the SZE/X-ray derived distances. The central values of the 
age of the universe derived from other methods as well differ from PDG value. For example 
to derived by the age determinations of elements by radioactive decay ratio method give the 



age of Milky Way ranging from 12.3 to 17.3 Gyr 12(|, the radioactive dating of old stars give 



values in the range 11 to 20.2 Gyrs 



13| . the age of the oldest star cluster ranges from 8.5 

to 16.3 Gyrs |l4|. Another point to mention is that the PDG value of to is determined from 

ACDM model. Therefore it is better to be more open minded to be about the value of Hoto, 

and hence in the following I take 

P2a2to 1 -i -i 

^0^0 = : + ^ha^ "tg = ~ 1 

Pi + P2a2to 6 

^ P2a2to = /3(e-^ln/3) , pi = /3(l-e + ^ln/3) (22) 

where Eqs. ([3]) and fl20|) are employed. 

One may obtain a constraint on the value of /3 by using the cosmic deceleration period 
(in the matter dominated era). Eq.(j5]) suggests that at the matter dominated era 

1 f h \ ^ 

-7piia2tm)^ + Pibi < t^> — fyj (23) 

-7 Pi + hp2a2t^ < ^ t^ < (24) 

OP2(l2 

1 [hi\' 7pi 



0.2 \ 7 J 5^20-2 

1 V 7(l-^ + iln/3) , , 

7'"^) <^"< skIim) '''' 

where tm denotes the time of deceleration in the matter dominated era, and 7m = 
Note that the inequalities above do not saturate i.e. the lower and the upper values in the 
inequalities are not infinitesimally close to the initial and final times of cosmic deceleration. 
In fact a more stringent bound on /3 and the time of the onset of cosmic acceleration in the 
dark energy dominated era may be obtained . It is evident from ^ that t' = is greater 
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than the time of onset of dark energy dominated era, i.e. = at' , a < 1 because of the 
additional terms contributing to the denominator of Eq. in addition to those considered 
in Eqs.f l23|l24|) . Hence at the onset of cosmic acceleration one may write 



bi{pi + P2a2td) + {ci2td) 6 (-7pi + 5p2a2td) = 

bi{Pi +P2a2C(t') + {a2at')^{-7pi + 5p2a2at') 

bi 5 1 

= y(7pi - 5^20-2^') + bi[-p2a2t' + P2a2t'] + {a2at')6[-7pi + 5p2a2t' + 5p2a2{a - l)t'] 

= bi{^ + a)p2a2t' + 5p2a2(a - l)t'{a2at')6 = 
=> + a) + 5(a - l)(a2at')^ = 

^ + «) + 5(a - 1)4 =0 ^ ln/3 = ^^^tI (26) 

where ■jd = ¥, and t' is the time satisfying —7pi + 5p2a2t' = 0. We know that td < to- 
The observational data analyzed in the context of ACDM model and dynamical dark energy 
models with a moderate dependence on redshift gives ^ | [15]. The fact that there is no 
significant disagreement of the ACDM model with data implies that the value of td should 



not be too different from this value. If one takes ^ = i(|^)^ 0.89 while for = 

to 2 \ to / to iUU 



1 



1 ,1 



^ ~ 0.464. Therefore it is safe to say that = (j^^ ^ ~ 1 for reasonable values of 7^. 
Then one may get an idea of the magnitude of In /3 for a few values of a by using Eq. fl26|) 



a = 1 =^ InP = 13 = 1 (27) 

/n/3 = ^7| => /3 ~ 1.63 (28) 



a = — ^ 

10 113 " 113 

a = ^ lnf3 = ^71 ~ ^ ^ /3 ~ 7.8 (29) 

44 



1 31 5 1 31 5 

a = - ln^ = ^7| ~ ^ => /5 ~ 251 (30) 

a = = I7I ~ 7 ^ /3 ~ 1097 (31) 

7 

It is evident that in any case 

ln/3 < 7 (32) 

In the following paragraphs we take this as an upper bound on the values of In /3 and we do 
not consider higher values unless it seems necessary for the sake of completeness. 

Now we derive a lower bound on the value of (3 by using the Goo/f at present time. 
Note that we use Goo/f rather than the equation of state for dark energy since the dark 
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energy and dark matter fields are mixed in the energy-momentum tensor so that it becomes 
impossible to entangle the dark energy and dark matter contributions properly in this case. 



' ^ |87rG(p + 3p)y bi{pi + p2a2to)[bi{pi + p2a2to) + {a2to)6 {~7p^ + 5p2a2to)] 



a ' t=to ^ o ^' t=to 

?2 



3 \ 36^' 



lnl3 J (ln/3 + 7-120 

The PDG values -1.14 < u^ark-e < -0.95, 0.21 < fi„ = < 0.26, Qdark-e = t^t - 
0.74 may be used to calculate f l55|) . The corresponding observational value of the ratio is 

3.9 < = < 5.5 (34) 

However there are studies with a wider range for current equation of state and density 
parameter from the analysis of SNe data alone |l^ 

-1.7 < codark-e < -0.50 , 0.23 < fi^ = — < 0.37 (35) 

Pm 

-6/(l + 3(-l,7)0.77) = 2,0B < (%)^^^^ = < oo (36) 

Although the infinity is unphysical I do not know a stringent and definite upper bound to 
be replaced by the oo in (1361) . Therefore I keep it as infinity. However one may replace oo 
by a large enough value. For example gives upper bound /3 ^ 23. We plot f %^') 

V a / t=to 

versus /3 for various values of 7 and ,^ ~ 1. The results are given in Table [H One sees that 
the values of /3 compatible with (IMI) are greater than 2.2 while the lower bound on /3 for 
fl5^ with (3u = 9 are greater than 2.6. 

As a complimentary analysis one may determine the ratio Gij/GoQ. Consider ^/Goo at 
present time 

^\ _ 12pi(a2to)i + bjjpi + p2a2tof + |&i(a2to)^ (-7p^ + 4pip2Q2to + Halpltl) 
^^00 J ^^^^ [6^2(02^0)^ + bi{pi +P2a2to)Y 

18^2_(7_12^)1^ ^^(1^^)2 ^^^^ 



54^ 



A (—/Goo) versus /3 graph may be plotted for various values of ^. I take —0.315 < 

V / t=to 

^dark-e^dark-e < —1.3629 by usiug Eq. (!3H|) . The corresponding allowed range of values of 



(33) 
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V " / t=tn 




0.5 


3.9 — 5.5 


none 


0.5 


2.05 — 23 


2.6 — 3 


0.8 


3.9 — 5.5 


none 


0.8 


2.05 — 23 


2.14 — 12.2 


0.9 


3.9 — 5.5 


none 


0.9 


2.05 — 23 


1.89 — 43.68 


1 


3.9 — 5.5 


3.95 — 144.4 


1 


2.05 — 23 


1.804 — 147.4 


1.2 


3.9 — 5.5 


3.55 — 1632.16 


1.2 


2.05 — 23 


1.757 — 1635.05 


1.5 


3.9 — 5.5 


3.56 — 4.34 


1.5 


2.05 ~ 23 


1.77 — 6.8 



TABLE I. The allowed values of (3 for two intervals of (Goo/f )^_^ ,and various values of ^ 



/3 for some values of ^ are given below 



e = 


0.8 


(3 


= 0. 


8 


16.8 


e = 


0.9 


/3 


= 0. 


82 


54 


e = 


1 


/3 = 


0.8 




-- 180 


e = 


1.2 


/3 


= 0. 


83 


2 X 10^ 



(38) 



In fact we should exclude the values of /3 smaller than one given above because of the 
definition of (3 in ( 120|) . The values of f3 above are barely consistent with the more stringent 
bounds in Table [T] for C, = 0.8, 0.9 and are consistent in the upper range for the others. The 
^^/Gooj value corresponds to the effective equation of state of the dark fluid consisting 
of dark energy and dark matter. Therefore it is useful to give its general time dependence 
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as well. 



FIG. 1. (f^/Goo) versus 7 = ^ graph for /? = 50, ^ = 1 



Gil 
911 



00 



[6^2(02^)6 + bi{Pi + P2a2tW 



127^ [l-(l-7)5+^ln/3]2 



+ 



ll((e-iln/3)2ln/ 



7(1-5+1 In /3) In /3 



187S(l-(l-5)7+ 



(1-7) 



"TT" 



ln/3) 187"B-(l-(l-5)7+ 



(1-7) 



ln/3) 



3(7-^¥ + 



5-1 In /3 



6 (l-(l-07+^ln/3) 



(39) 



where 7 = ^- As we shall remark later in this section a general analysis of this effective 
equation state at an arbitrary redshift is quite difficult due to the highly nonlinear form of 
the above equation. However one get an idea of its general variation by the inspection of 
Figured for /3 = 50, ^ = 1. 

The general form of the cosmic history must have a cosmic acceleration era corresponding 
to the time of inflation that is followed by an era of deceleration at the matter dominated 
era, and finally by the present time acceleration era. Moreover the redshift values and ages 
for these eras must coincide with the observational data jll] at least at the order of magnitude 
to have at least an approximately realistic model. For this purpose we draw f /Cqo versus 
time for C=0.8„ 1, 1.2; (3=2, 5, 10, 20, 50, 100, 200, 500, 1000, 2000, 10000, 3000, 40000 by 
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using 



FIG. 2. I/Goo versus 7 = f graph for /3 = 50, ^ = 1 



(02)661(^1 +^2^2^) [7-^ (Pi +P2a2t) + (-7^1 + 5^2^2^)] 



Goo 3{a2t) 3 [6^202^ H ^(Pi + ^20.2^)] 

{a2t)T> 

InP \ A[ij-l In /3) A - 7(1 - e + I In /3) + 57^ - §7 In P] 



37^ ; 6(^ - I In /3)7 + 7-hn /3[1 - (1 - 7)e + V 1^ /^l 



(40) 



A = [l-(l-7)e + ^l^ln/3] 

where j = j^. One may get a sense of the general form of the evolution in Figure [2] for 
/3 = 50, ^ = 1. I do not give the other plots not to make the paper too crowded. The 
all values of the times of the start and the end of cosmic deceleration that are not wholly 
excluded by data are given in Table HTl below. It seems that this analysis prefers lower values 
of p. Note that we should not expect a good match between the values obtained here and 
the observational data for the time of the start of the deceleration period because at that 
time the radiation has a major contribution and we neglect the contribution of radiation in 
this study. 

Now we compare the data and the predictions of this model for different redshifts and 
times, and hope, at least, an order of magnitude agreement. First consider the time of 
the starting of cosmic acceleration. This is the time where the deceleration changes into 
acceleration, hence the acceleration of the cosmic expansion is zero i.e. the numerator of ([5]) 
is zero, namely 

/ = + [1 - (1 - 7.)e - 14 - llh^ - 7(1 - O7I = (41) 



Id = ^ , X = InP 
to 



(42) 
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/3 








P 






0.8 


2 


10-6 


0.5 


0.8 


5 


1.5 X lO-'' 


0.7 


0.8 


10 


1.2 X 10-3 


0.9 


0.8 


200 


0.165 


0.947 


0.8 


500 


0.28 


0.89 


0.8 


1000 


0.35 


0.867 


0.8 


2000 


0.58 


0.91 










1 


2 


9.5 X lO-'^ 


0.13 


1 


5 


1.4 X 10-^ 


0.26 


1 


10 


1.3 X 10-3 


0.34 


1 


20 


6.6 X 10-3 


0.415 


1 


50 


0.04 


0.45 


1 


100 


0.116 


0.7 


1 


1000 


0.48 


0.945 


1 


2000 


0.52 


0.916 



TABLE II. The times of start, ^-f^ and end, of the cosmic deceleration for different values of /? 

Because of its highly non-hnear form this equation I could not analytically solve this equa- 
tion. However after plotting / versus x for various values of ^ and 7^ and determining the 
location of zeros one may get some information. The result is given in Table IIIII Keeping 





Id 


X = ln/3 




Id 


X = ln/3 


0.8 


0.85 


2 or 26 


0.8 


0.8 


2 or 25 


0.8 


0.75 


2 or 20 


0.8 


0.7 


2 or 16 


0.8 


0.5 


0.7 or 9 


0.8 


0.4 


0.2 or 7.8 


0.8 


0.1 


1 or 5 


0.8 


10-10 


-1.2 or 0.2 


0.8 


10-20 


-1.2 or 


1 


0.85 


5 or 35 


1 


0.8 


5 or 24 


1 


0.75 


5 or 17 


1 


0.7 


5 or 15 


1 


0.5 


3.6 or 7.2 


1 


0.3 


1.9 or 5.5 


1 


0.1 


0.5 or 4.5 


1 


10-10 


0.06 or 0.09 


1 


10-20 


or 0.0033 


1.2 


0.85 


7 or 33 


1.2 


0.8 


7 or 22 


1.2 


0.75 


8.7 or 15.7 


1.2 


0.5 


none 



TABLE III. The allowed values of lnj3 for various values of ^ and "fd 

these values in mind now we may find the redshift values and the time of onset of current 
cosmic acceleration, td predicted by this model and compare the observational values given 
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in literature. Consider aQ/a{tf}) 

Qq ^ exp [&i(a2trf)~g] ^ /g'^^^" ^ .^g. 

a(td) (Pi+P2a2td) (l-0 + |ln/3 + (e-|ln/3)7d ^ ^ 

The analysis of cosmic data [18] gives the redshift and time of onset of dark energy dominated 
era, respectively, in the ranges z = 0.66 — 1.21. td = (5.7 — 8.5) Gyr. The allowed intervals 

of P in fj43|) where is in the range 1.66 2.21 for the phenomenologically relevant 

values of ^ and 7d = 7^ for ^= 0.8, 1, 1.2 and 7^=0.1, 0.4, 0.5, 0.6, 0.7, 0.8 may be found in 
Table HVl Table tells us that 7^ = 0.1 is inconsistent with data. Comparison with Eg. (138 p 





Id 


/? 


0.8 


0.1 


none 


0.8 


0.4 


0.15 — 123 


0.8 


0.5 


1 — 7 X 10^ 


0.8 


0.6 


2 X 10^ 2.1 X 10^ and 1.8 x 10"^ 2.5 x 10"^ 


0.8 


0.8 


5 X 10^2 2.5 X 10^^ and 7 x 10^^ 3 x lO"'^ 


1 


0.1 


none 


1 


0.4 


none 


1 


0.5 


4 — 1400 


1 


0.6 


130 7 X 10^ and 1.2 x lO^^ l 


1 


0.7 


7.5 X 10^ 10^° and 1.7 x 10"^ 5 x lO^^ 


1 


0.8 


1.5 X 10^2 iQi8 j^j^j 5 X 2.8 X 10-6 


1.2 


0.1 


none 


1.2 


0.4 


none 


1.2 


0.5 


barely 30 


1.2 


0.6 


0.5 2 X 10^ 


1.2 


0.7 


9 X 10^ 3.5 X 10^ and lO'^ 9 x lO'^ 


1.2 


0.8 


10^2 5 X 10^7 j^j^jj 2.6 X lO"'^ 2.1 X 10-5 



TABLE IV. The allowed range of values of (3 for various values of ^, 7 with -fP-r in the range 
1.66 - -2.21 

and Table H] implies that the 7^ values in the range 0.4 - 0.6 are consistent with redshift data 
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18| and the time of the onset of the cosmic acceleration for the phenomenologically relevant 
values of in the range 1.2 - 0.8. 

An important point is to be mentioned at this point: Note that the values 2; = 0.66 — 1.21, 



td = (5.7 — 8.5) Gyr in 15j are derived by the assumption that the Hubble constant at 



scale factor a{t) may be expressed as 

H{a) = Ho[—^ + ^ + ^+ f.. 2 44) 

In principle one may define an effective equation of state u^ff as in [19] when dark matter 
and dark energy are coupled. However in this model the contributions of dark matter and 
dark energy are not only coupled they are mixed. Therefore their contributions can not be 
separated from each other properly. Moreover dark matter in this model is not dust-like (it 
only mimics a dust in the matter dominated era) while the matter in the above equation 
is dust-like. Furthermore in 15| and similar studies an equation of state for dark energy of 
the form u}{a) = —uq + Ui{l — a) or similar forms are employed. Let alone that a proper 
equation of state for dark energy in this model can not be defined a common equation of 
state for dark energy and dark matter is highly nonlinear as seen before in Eq. fl39p 



^ _ [l-(1^7)g+V^ln/3]^ i8^^(l_(l_g)^+(l^ln/3) ig^^ ln/3) 

' + (l-(l-07+V^ln/ 

An inspection of [a(t)]^^ versus 7o = ^ graphs show that in the low redshift range z= 05 
- 2 one may approximately take [a(t)]^^ proportional to 7. In other words one may get 
the form of uj{a) for low redshifts by simply replacing 7 in fHSjl by 1+z. It is evident that 
this relation is quite nonlinear in z. For higher redshift values the relation between 7 and 
(z+l)~^ also becomes non-linear making the form of u{a) even more complicated. Therefore 

;his model with data in a more precise way 
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18| with keeping these points in mind. 



in order to see the degree of the compatibility of ; 
it is necessary to repeat the analysis of data in 
Only then one can say some definite conclusion on the degree of the agreement between this 
model and observational data. In any case I think the rough analysis given in this paper 
is enough to consider this model as viable toy model in the direction of unification of all 
eras of cosmic history. In fact all I have mentioned in the context of the analysis of the 
onset of cosmic acceleration data is true for the analysis of data for equation of state of dark 
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energy 



20| . density parameters |16[. and the analysis of data on time reionization and time 
of matter-radiation decoupling |21[ discussed below. 

Before continuing the comparison of the model with observational analysis for the times 
of reionization and decoupling now I want to consider the inflationary era because there is 
no baryonic matter or radiation effect in this era, this toy model is expected to be most 
similar to the reality in this era in the context of this model. It is evident from (jS]) that, at 
the time of inflation, 

biiPi +P2a2t) + {a2t)^{-lpi + 5^202^) > (46) 

This condition is satisfied for very large and very small a2t's. We identify the very small 
values that satisfy fH6|) as inflationary times, and at small times fH6|l is guaranteed if we 
take a2ti <^ 1 where the subindex i refers to inflation. In fact a more stringent bound may 
be obtained from the slow-roll parameter in (|4]) 

where we have used the fact that the pi term in (jlj) is the leading term in the inflationary 
period. Then 

ajtis) _ jpi + P2a2Us) exp {-bi[{a2tisy^ - {a2tiiy^} 
aiUi) {Pi+P2a2tii 



^. (1-7.), n I exp{- — [^./ -hi']} ^ exp{- — [t,/ -t,1} 

(48) 

where t^s and tn are the times of the start and end of inflation, respectively. If we assume 
60 e-fold expansion and t^g = = 10~^^ sec, tn = 10^^^ sec 

^ exp [-61 {a2)--^lQ\sec)--^] = e*^" 



a{tii 

=^ bia2^ X lOSec-H = 60 

=^ 6ia^^ ~ 6x lO^^sec^ ln/3 = 8.8 x 10"^ (49) 

If we assume 60 e-fold expansion and tis = = 10"^*^ sec, Ui = tig + lO"^"^ sec 
"'^''^ - exp {-61 (02)^^(1 X 10"^°)"^[(100)"^ _ (100.0001)"^]} = 



=^ 6ia^^ ~ 7.8 X 10^ ^ ln/3 = 11 (50) 



1 
1 



feitto " X 7.7 X 10"^ = 60 
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Note that 

6102^ = bi{a2to)~hl = 4 ln/3 (51) 
Some other values of t^s, tu and ln/3 for 60 e-fold expansion are 

Us = 10"^° sec Ui = tis + lQ-^'^ ln/3 = 5.57 x 10"^ (52) 

= IQ-^Sec Ui = ti, + 10-^Sec ln/3 = 5.57 x 10^ (53) 

t.^ = 10-28 sec + IQ-^Sec tu = t,, + 2 x lO^^Sec ln/3 = 5.72 (54) 
Us = 10^2^ sec + IQ-^Sec Ui = Us + 2A x 10~^Sec ln/3 = 8.15 (55) 

t.^ = 10-29 sec Ui = Us + 5 X 10-^^ sec In /3 = 0.13 (56) 



One notices that (H9|) . ( 152|) . ( I54|) . (156|) are consistent with ( l32|) while the others are not. 
However it seems that the values in Table I exclude the values of In /3 much smaller than 1. 
This excludes the options in (H^ and as well. Hence the viable values seem to be 
and (jSSD and all values of parameters between them and close to these values. This offers a 
wide range of tis between IQ-^* sec and IQ-^^ sec. It is evident that all phenomenologically 
viable values may be obtained by adjusting Us the tis = 10~^^ sec - tis = lO^^s _gg^ range 
that corresponds to a lower scale inflation ^22] . 

A comment is in order at this point. From Eq.© we see that just at the end of the 
inflationary era 

bi{Pi +P2a2tmi) + {a2tmi)H-7Pi + ^P2a2tmi) < (57) 
Pi— ^-7pi + e = (58) 

[0,2tmi)^ 

(7J ln/3 - 7) (1 - e + ^ ln/3) + e' = (59) 
—P2a2t + 5p2a2t , e' = = ImiilrJ + 5)(^ - ^ ln/3) 



{a2t)6 P 6 

where the first two terms in f l58p are the dominant terms and e (and e') is small with respect 
to the others. The fact that e' in and at the end of inflationary era is small implies that 

either ln/3 - 7) or (1 - ^ + | ln/3) is small. Taking tmi ~ tu ~ IQ-^^ IQ-^^ i.e. 

_i 

Imi = ^ ~ IQ-^^ implies that (7^/ In /3 — 7) is not small unless ln/3 is extremely small. 
Therefore (1 — ^ + | In /3) should be small if deceleration era starts just after the inflationary 
era. This may be provided by taking ^ a little bit larger than 1 and In /3 small. For example 
one may take ^ = 1.05 and ln/3 ~ 0.3 (i.e. (3 ~ 1.35). Otherwise one should take the start 
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of the deceleration era much later than the standard inflationary era (i.e. the inflationary 
era is much longer than the standard inflationary times). Although this option seems to be 
a less acceptable option it is, in fact, the more reasonable choice. This is due to the fact that 
we neglect radiation in this study. In the realistic case there is a radiation dominated era just 
after the inflationary era. Radiation like matter drives the universe towards deceleration. 
Therefore if we add radiation to the model it is effect will be an earlier start of deceleration 
era compared to the radiationless case. This explains why the time of the start of the 
deceleration period almost coincides with the time of start of the matter dominated era in 
Table unless (3 is extremely close to 1. In other words the values of parameters become 
less reliable as we get closer to the radiation dominated era. We should keep this in mind 
as we analyze the observational data. 

Now we apply the values obtained to the time of reionization, t^i- In fact we expect, 
at most, a rough agreement with data since tri goes deeper into the matter dominated era 

where neglecting baryonic matter becomes more questionable. 

^ _ 1 

ao exp [6i(a2tri)"«] Z?"^" 



aiUi) iPi+P2a2tri) (l-0/3 + |ln/3+(e/3-|ln/3)7. 



(60) 



where jri = j^- The observational value of is 12 ± 1.4 and the corresponding ACDM 

value of tri is 430^70 that corresponds to jri in the interval 0.0225 0.0472 if 

one assumes a loose bound on the value of to; = H 16 Gyrs in the light of 

the values of to from different observations mentioned before. One may plot versus 
7 = ^ for the phenomenologically relevant values of ^ and various (3 values. The allowed 
intervals of 7^ for ^ in the interval 10.6 - 12.4 for various values of ^ and /3 may be 
found in Table |Vl It seems that the values of (3 compatible with data are 10, 20 for ^= 
0.8; 2, 5, 10 for ^= 1; and none for ^=1.2. However one should keep in mind that a more 
detailed analysis may give a wider range of parameters since the age calculations in the data 
analysis [l3| use a restricted form for dark energy where Hubble constant may be expressed 
in terms of density parameters where matter is assumed dust-like and a restricted form of 
variation of dark energy with redshift, and a restricted class of equations of state for dark 
energy where dark energy is not entangled with matter as pointed out before. Therefore 
reanalysis of data in the context of this model is necessary to reach a more precise and more 
definite conclusion. Another factor for poorer agreement with data is that we neglect the 
contribution of baryonic matter whose contribution in matter dominated period is greater. 
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/? 


7ri 




/3 


7ri 


0.8 


2 


0.0016 ^ 0.0023 


0.8 


5 


0.0115 ~ 0.016 


0.8 


10 


0.023 ~ 0.03 


0.8 


20 


0.037 — 0.046 


0.8 


50 


0.058 — 0.068 


0.8 


100 


0.074 — 0.086 


0.8 


200 


0.091 — 0.105 


0.8 


500 


0.116 — 0.128 


0.8 


1000 


0.13 — 0.146 


0.8 


2000 


0.15 — 0.162 


1 


2 


0.028 — 0.041 


1 


5 


0.029 — 0.041 


1 


10 


0.04 — 0.051 


1 


20 


0.054 — 0.066 


1 


50 


0.079 — 0.086 


1 


100 


0.089 — 0.104 


1 


200 


0.105 — 0.12 


1 


500 


0.13 — 0.142 


1 


1000 


0.142 — 0.16 


1 


2000 


0.16 — 0.18 


1.2 


2 


0.17 — 0.19 


1.2 


5 


0.105 — 0.125 


1.2 


10 


0.09 — 0.105 


1.2 


20 


0.091 —0.104 


1.2 


50 


0.098 — 0.116 


1.2 


100 


0.114 — 0.127 


1.2 


200 


0.125 — 0.142 


1.2 


500 


0.144 — 0.162 


1.2 


1000 


0.159 — 0.177 


1.2 


2000 


0.173 — 0.192 



TABLE V. The ahowed range of values of "yri = ^ for various values of ^, 7 with ^^^^ in the range 
10.6 - - 12.4 



Next consider the data for the time of decoupling and the corresponding redshift; ~ 
1090 

_ 1 

ao _ exp [6i(a2t*)"^] _ /3^* ^ ^^^^ 



where 7^, = ^. One may plot versus 7 = ^ foi' various values of C, and /3. The 
values of 7=,. corresponding to the observational value of ~ 1090 are given in Table 
IVII Inspection of the table suggests there are no of the values are ^ and /3 compatible 



with observational value 3.8 x 10^ yr 17| that corresponds to the interval 7* = 



2.375 X 10"^ 3.4545 x 10^^ (provided that to = (H 16) Gyrs) except for 

^ = 0.8, 1 and /3 somewhere between 2 and 5 (i.e. ~ 3.7, 3) while the values for ^ = 0.8, 
/3 = 5 and ^ = 1.2, /3 = 5 are close to the relevant values. In fact this poor agreement with 
those given in [itI is expected. In addition to the reasons mentioned for the reionization 
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Iri 




/3 


Iri 


0.8 


2 


1.42 X lO"'^ 


0.8 


5 


7.45 X 10-^ 


0.8 


10 


3.3 X 10--^ 


0.8 


20 


9.05 X 10-"^ 


0.8 


50 


2.32 X 10-3 


0.8 


100 


3.99 X 10-3 


0.8 


200 


6.17 X 10-3 


0.8 


500 


9.86 X 10-3 


0.8 


1000 


1.32 X 10"2 


0.8 


2000 


1.71 X 10-2 


1 


2 


3.9 X 10-6 


1 


5 


1.15 X 10-"^ 


1 


10 


4.44 X 10-"^ 


1 


20 


1.12 X 10-3 


1 


50 


2.7 X 10-3 


1 


100 


4.5 X 10-3 


1 


200 


6.8 X 10-3 


1 


500 


1.06 X 10-2 


1 


1000 


1.425 X 10-2 


1 


2000 


1.81 X 10-2 


1.2 


2 


7.91 X 10-2 


1.2 


5 


7.45 X 10-^ 


1.2 


10 


3.3 X 10"^ 


1.2 


20 


1.56 X 10-3 


1.2 


50 


3.34 X 10-3 


1.2 


100 


5.28 X 10-3 


1.2 


200 


7.7 X 10-3 


1.2 


500 


1.17 X 10-2 


1.2 


1000 


1.53 X 10-2 


1.2 


2000 


1.95 X 10-2 



TABLE VI. The allowed range of values of 7^, = for various values of ^, 7 with ^^^^ ~ 1090 

time there is an important additional source of discrepancy. The time of matter radiation 

decouphng time is quite close to the radiation dominated era. The ratio of radiation in this 

period in the order of a fourth of the total energy density at this time while this model 

neglects the contribution of radiation. 

To summarize the results of this section can be stated as follows: We have seen that 

the predictions of this model for each of (-/Goo)._. , (—/Coo) ? o-o/o-itd)-, CLo/o-iUi) 

« t-to \ gii y j^^^ 

are compatible with observations although not with central values given in literature. The 
prediction of the model for ao/a(t*) is partially consistent with observational values. In fact 
the relatively less compatibility for the decoupling time is expected since the radiation- 
matter decoupling time is close to the radiation dominated era while radiation is ignored 
in this study. I have also shown that an inflationary era naturally fits the model. One 
may consider the simultaneous compatibility of the predictions of all these parameters with 
observations as well. In all cases there is wide range of /3's compatible with Eg. (15^ . The 
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values of /3 allowed by (1551) includes the values allowed by Table HI that is, Table H] and 
f l55]) are compatible while Table H] is more restrictive. The values of 7^ in Table IIIII that are 
compatible with Table HV] are ^=0.8 7^=0.4 or 7^=0.5, ^=1 7d=0.5 or 7^=0.8, ^=1.2 
=^ hardly 7^=0.8. The values of ^, 7^ in Table IIVI whose (3 values compatible with the /3 
values in Table E are ^=0.8, 7d=0.4, 7^=0.5, (3= 2.14 - 12.2; ^=1, 7d=0.5 ^ /3= 4 - 144; 
^=1, 7^=0.6 =^ f3= 130 - 147; ^=1.2, 7^=0.6 =^ Table|I]and Table HVl are compatible except 
for lowest values of /3. Mostly Table [T] is more restrictive than Table HVl The values of 'jri in 
Table rVl that are more compatible with observational value 7ri= 0.0225 - 0.0472 in literature 



171 ] seem to prefer /3 in the range 2-10. I do not use Table IVT] to constraint (3 since the time 



of decoupling is close to the radiation dominated era while we do ignore radiation, so the 
reliability of the values obtained is questionable, and an order of magnitude compatibility 
is enough. We see that compatibility of the values of all these tables seem to prefer values 
^ = 0.8 - 1, /3 = 2 - 10. However these values of f3 are at the edge of the observationally 
allowed values rather than being centrally allowed values. The limited overall compatibility 
of the results of this model with observations may either be due to this model being simply 
a toy model or the inapplicability of some of the assumptions of the analysis in literature 
to this model such as Eq. fj4^ and uj{a) = —Uq + (^i(1 — a) or a combination of both. In 
fact, even a standard analysis may be enough to check the viability of this model beyond a 
toy model for small enough redshift bins. For example, it seems that the allowed value o" 
the equation of state for dark energy, ude the smallest redshift bin in Figure 14 of {20 1 
may be as large as -1/3 while the Gn/Goo versus f3 graph for this model at present time 
(for ^ = 1) gives Gn/Goo — —0.45 (that corresponds to ude ~ —0.65 for Qde = 0.74) for 
/3 ~ 3. A definite conclusion needs a detailed comprehensive reanalysis of all data in the 
light of this model in a separate study. 



IV. CONCLUSION 

I have considered a model where inflationary era and (dark) matter dominated eras are 
induced by a scalar field 0i while the dark energy dominated era is induced by another scalar 
02- These fields may be either considered to be fundamental fields or as effective classical 
fields. I prefer to consider them as classical fields rather than true fundamental fields. I have 
neglected the effects of baryonic matter and radiation. A rough phenomenological analysis 
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of cosmic data gives an order of magnitude agreement with data. This is encouraging for 
future studies in this direction. One must include baryonic matter and radiation to obtain 
a more reahstic model. However this is not an easy task. First difficulty is that baryonic 
matter and radiation should be included after the time of inflation because it should be 
produced by the decay of one of the scalars (probably by (pi). Second, even when one 
includes them in ad hoc way this modifles the metric. Hence one must flnd the scale factor 
that corresponds to inclusion of the baryonic matter and radiation and this not a trivial 
task. Another point that needs further study is a more detailed and comprehensive analysis 
of the available parameter space and to flnd the most optimal set. Yet another point for 
further study is the study of cosmological perturbations produced in the inflationary epoch. 
The inflation obtained here is a standard slow-roll inflation with the canonical kinetic terms 
for the scalars. Therefore the general form of the perturbations is the same as the usual 
slow-roll case 23|. However a detailed study of the perturbations in this model should be 



obtained to compare with the expectations of the other models for data to be obtained in 
future cosmological observations. All these points need further separate studies. 
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